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ABSTRACT 



Context. An increasing number of potentially habitable terrestrial planets and planet candidates are found by ongoing planet search 
programs. The search for atmospheric signatures to establish planetary habitability and the presence of life might be possible in the 
future. 

Aims. We want to quantify the accuracy of retrieved atmospheric parameters (composition, temperature, pressure) which might be 
obtained from infrared emission spectroscopy. 

Methods. We use synthetic observations of the atmospheres of hypothetical potentially habitable planets. These were constructed with 
a parametrized atmosphere model, a high-resolution radiative transfer model and a simplified noise model. The simulated observations 
were used to fit the model parameters. Furthermore, classic statistical tools such as x 2 statistics and least-square fits were used to 
analyze the simulated observations. 

Results. When adopting the design of currently planned or proposed exoplanet characterization missions, we find that emission 
spectroscopy could provide weak limits on surface conditions of terrestrial planets, hence their potential habitability. However, these 
mission designs are unlikely to allow to characterize the composition of the atmosphere of a habitable planet, even though CO2 is 
detected. Upon increasing the signal-to-noise ratios by about a factor of 2-5 (depending on spectral resolution) compared to current 
mission designs, the CO2 content could be characterized to within two orders of magnitude. The detection of the O3 biosignature 
remains marginal. The atmospheric temperature structure could not be constrained. Therefore, a full atmospheric characterization 
seems to be beyond the capabilities of such missions when using only emission spectroscopy during secondary eclipse or target visits. 
Other methods such as transmission spectroscopy or orbital photometry are probably needed in order to give additional constraints 
and break degeneracies. 

Key words. Planets and satellites: atmospheres, techniques: spectroscopic, methods: data analysis 



1. Introduction 

The possibility of finding potentially habitable or even inhab- 
ited terrestrial planets is one of the exciting motivations for the 
search of extrasolar planets. So far, more than 50 planets with 
(minimum) masses below ten Earth masses are known among 
the more than 800 detected extrasolar planets. Furthermore, re- 
cent studies based on radial-velocity planet searches claimed 
that the occurrenc e of low-mass, potentially terrestrial planets i s 
rather large (e.g., iHoward et aDl2010l IWittenmver et al.ll201 ll) . 
For short-period orbits with periods < 50 days, 11-17% of all 
stars host at least one low-mass planet. Estimates inferred from 
microlensing surveys i mply that the mean number of planets per 
star is larger than one (ICassan et al.ll2012T) . 

Already, some potentially habitable (candidate) super-Earths 
in or very close to the habit able zone (HZ) of their cen- 
tral s t ar have been discovered (lUdry et al.l 2007[ Mayor et all 
20091 iBorucki et al.ll201 il iPepe et al.ll201 ll iBonfils et alJl2012l 



biomarkers are necessarily surface or atmospheric spectral sig- 
natures. 

Several studies h ave addressed the a spect of atmospheric 
biosignatures (e.g., Sagan et al.l 1993L Schindler & Kastingl 



1200(1 iDes Marais et al.l 120021 ISelsis et al.l 120021 ) or the re 
sponse of spectra with respect to central star and atmo- 
spheric composition (e.g., Segura et aD 20031 ISegura et al.l2005l 



iKaltenegger et al.l2007llVasquez et al.l2013 ). Further studies fo 



tectability (e.g., 


Kaltenegger & Traubl2009l Deming et al.ll2009L 


lBeluetalJl2011 


.iRauer et al.ll201 ll von Paris et al.ll201 lb. These 



Anglada-Escude et al.ll2qi2LlDelfosse et al.[l2oT2l IBorucki etail 
20121 iTuomi etalJ2013 ). Hence, the detection of potentially 
habitable terrestrial planets is within reach of present-day in- 
strumentation. The next step would be to search such potentially 
habitable worlds for the indications of the presence of life, so- 
called biomarkers. To be detectable remotely from Earth, such 



studies calculated signal-to-noise ratios and suggested observa- 
tion strategies using instrument capabilities of the planned James 
Webb Space Telescope (JWST). Since exoplanet characteriza- 
tion is not the design-driving purpose of JWST, dedicated ex- 
oplanet space mission concepts have been developed. Some of 
them have already been proposed to ESA or NASA. The con- 
cepts include coronographs such as SEE-COAS T or ACCESS 
(e.g., ISchneider et all l2009l iTrauger et al 1 120081). interferome- 



ters such as Darwin ( e.g. jLeger et al. H996tlCockell et al.l2009al 
ICockell et al.l l2009bh and ne ar- to mid-IR spectrographs such 
as EChO or FINESSE (e.g., iTinetti et al.l I2Q12L iTessenvi et all 
l2^ISwair3l2012h . The ultimate aim of these missions is to 
characterize the atmospheres of (habitable) exoplanets spectro- 
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Table 1. Emission spectroscopy of potentially habitable planets: Science goals (species to be detected and characterized) and 
requirements in terms of spectral resolution and S/N ratio for two European exoplane t characterization missions. For Darwin, these 
S/N goals limit the number of potential targets to 40 (20) for R=5 (R=20) (see, e.g.. lLeger et al.lfl996b . In the case of EChO, tens 
to hundreds of tran sits need to be accumulated in order to achieve the stated aims, depending on the IR magnitude of the host star 
dTinettiet al.ll2012l) . Because of mission lifetime constraints, EChO is limited to M-type stars brighter than about K^9 m . Currently, 
there are only three transit host stars satisfying these criteria, all orbited by hot Neptunes or hot super-Earths. 



Mission 


Science requirement 


Goal resolution 


Goal S/N ratio 


References 


EChO 


H 2 0, C0 2 , 0, 


10 


5 


Tinetti et al. (2012),Tessenyi et al. (2012) 


Darwin 


C0 2 , 6 3 


20 


5 


Cockell et al. (20093) 


Darwin 


H z O, C0 2 , Ch 


20 


10 


Cockell et al. (2009b) 


Darwin 


C0 2 


5 


12.5 


Leeeretal. (1996) 


Darwin 


3 


20 


20 


Leger et al. (1996) 



scopically. The Darwin mission, which is not pursued further 
by ESA, was explicitly designed for the investigation of terres- 
trial, potentially habitable planets. The main focus of EChO, cur- 
rently in competition for the ESA M3 launch slot, is the obser- 
vation of hot Jupiter and hot Neptune planets. However, an addi- 
tional science opportunity is the possibility of observing (in case 
they would be discovere d) nearby (habitab le) terrestrial planets 
around low-mass stars dTinetti et al.l 120121) . Table [TJ states the 
current mission design of Darwin and EChO with the aimed 
spectral resolution and signal-to-noise (S/N) ratio for reaching 
the scientific goals and opportunities regarding terrestrial plan- 
ets. 

However, the difficulty of characterizing the atmospheres 
of small, potentially low-mass exoplanets even with high spec- 
tral resolution and good S/N ratios is illustrated by the case 
of GJ 1214b. Numerous studies using both ground- and space- 
based spectroscopy and spectrophotometry (e.g, Bean et al.l 



based spectroscopy and spectrophotometry (e.g, tsean et al 
2010l.lD6sert et al.ll201 lllCroll et al.l201 lUCrossfield et al.ll201 1 



Bean et al.ll20TTl Ide Mooii et al.ll2012llBerta et alj|2012t) aimed 
at constraining the range of possible atmospheric scenarios. 
Most observations favor a cloud-f ree water vapor atmo sphere. 
However, a data point presented by Ide Mooii et al.l d2012l) seems 
to favor a hydrogen-rich atmosphere with a possible low-altitude 
cloud layer. 

In view of such difficulties, several methods have been pro- 
posed recently to allow for an efficient search of sets of at- 
mospheric parameter s which produce a reasonable fit to the 
data. S tu dies by, e.g [Madhusudhan & Seagel (120091) . iLee et alj 
(120 121) or Line et al 1 (120121) focused on hot Jupiter planets, with 
both pri mary and secondary eclips e data from real observations, 
whereas iBenneke & Seager|(l2012l) investigated in detail the po- 
tential of transmission spectroscopy to characterize the atmo- 
spheres of hypothetical super-Earths. These studies have shown 
that it is in principle feasible to constrain atmospheric composi- 
tion and temperature structure for hot Jupiters and super-Earths, 
however given the current data quality, both in terms of spec- 
tral coverage and signal-to-noise ratios, error bars on retrieved 
atmospheric parameters are rather large. 

Such retrieval studies could provide valuable clues to the de- 
sign of the aforementioned space missions and should be taken 
into account when calculating integration times or allocating ob- 
servations. Therefore, we present here for the first time a study 
of the possible retrieval of atmospheric and planetary parame- 
ters for hypothetical habitable terrestrial planets from emission 
spectra. The focus of the study will be on possible constraints on 
surface conditions and detection of important molecules in the 
context of life, such as oz one and water. In contrast to t he work 
by, e.g.. ILee et al. I d2012l) or lBenneke&Seagerl (1201 2b . we do 
not use here a Bayesian optimal estimation approach but rather 



focus on straightforward \ 2 calculations and nonlinear least- 
squares fits to illustrate the difficulties of atmospheric retrieval 
for habitable planets. From there, we investigate quantitatively 
the potential of planned or near-future space instrumentation to 
characterize habitable planets. 

The paper is organized as follows: Section [2]presents atmo- 
spheric scenarios and models used. Results will be shown in 
Sect. [3] and discussed in Sect. [4] Conclusions are given in Sect. 



2. Models and atmospheric scenarios 

2.1. Parametric atmosphere model 

We use a semi-empirical parametric atmosphere model to de- 
scribe the temperature structure and chemical composition 
of an arbitrary terrestrial, cloud-free habitable atmosphere, 
simila r t o the appr o ach p resent ed in iMadhusudhan & Seagerl 
d2009l) . iLine et alj d2012l) or ILee et al l d2012l) . Note that 
simple, analytical ID mo dels have been proposed recently 
dRobinson & Catling! I2012I) which might also be applicable to 
the retrieval problem. 

The model atmosphere is divided into 50 layers, with a 
higher resolution in \og{p) in the lower atmosphere where the 
temperature gradient is expected to be larger. We make the fol- 
lowing assumptions regarding our temperature structure: 

- The lower atmosphere is convective, i.e. a troposphere ex- 
ists. The temperature gradient is given by the wet adiabatic 
lapse rate, taking into account the release of latent heat by 
condensation of water: 



d\n p (1 + jfY-cn 2 o)(l + /c H2 o) 

TkTT ~~ rwet ~~ rdf y ; TTTTi 



dlnT 



, fM lr L 2 

1 + TJu^^-O 



(1) 



where Fd iy = y~ (c p heat capacity of the atmosphere, R g!is 
the universal gas constant) is the dry adiabatic lapse rate in 
pressure coordinates, ch 2 o is the concentration of water, M w 
the molar mass of water, L the latent heat and f = M- the 
mass ratio between dry air and water. Throughout this work, 
we tookL=583cal g -1 , i.e. its value at 300 K. 
- Above the tropopause, the atmosphere is in radiative equi- 
librium. The temperature gradient is assumed to be linear in 
In p, given by a linear interpolation between the tropopause 
temperature and the temperature at the model lid. 

Four parameters are used for the temperature-pressure pro- 
file T = T(p). These are two temperatures, 7toa (temper- 
ature at the top of the atmosphere, TOA) and r sur f (surface 
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temperature) as well as two pressures, /? trop (location of the 
tropopause) and p^ (dry surface pressure). The total surface 
pressure p sm i is constructed from and the water vapor pres- 
sure /?vap,H 2 o(7 , surf) at the specified surface temperature r sur f, i.e. 
Psurf(Tsurf) = Pdry + Pva P ,H 2 o(r sll rf)- The TOA pressure p T0A 
is fixed at /?xoA=10~ 4 bar since the low- to medium-resolution 
emission spectra presented in this work are not likely to be sen- 
sitive to pressures much lower than that. Fig. [1] illustrates our 
parametric temperature-pressure profile. 

Parameterized temperature-pressure profile 





*>T /p 

J** TOA' TOA 






Stratosphere 




. < A p / Pirop 








Troposphere 








« T ,u/P 



T [K] 



Fig. 1. Illustration of the parameterized temperature structure. 

For ptmp^P^Psmi, the temperature r ( "' in a layer n (starting 
at the layer above the surface, with rc sur f ace =50) is calculated via 
(see eq.[TJ 



n(n+l) 



exp 



In 



(5?) 



(2) 



For pTOA<P<Ptrop, the temperature profile is calculated with 



T (n) = T T0A + 



Tu-op - Tjo\ 

V Ptoa I 



In 



Ptoa 



(3) 



where r trop (the temperature of the tropopause at p t rop) is calcu- 
lated based on the values of T SUI f, p mi f and p tTO p (see eq.[2j. 

The atmosphere is assumed to be composed of four species, 
i.e. CO2, H2O, the so-called biomarker O3 as well as N2. We 
chose CO2 and H2O because they are the main greenhouse gases 
on Venus, Earth and Mars, and both are generally accepted con- 
stituents of the atmospheres of potentially habitable planets. N2 
is a major atmospheric species for all terrestrial-like atmospheres 
in the solar system (i.e., Venus, Earth, Mars, Titan). We chose O3 
as the biomarker in our model atmospheres because it is usually 
considered to be a very prominent bioindicator (e.g.JSagan et aO 
119931 ISchindler & Kastingll2000l ISeTsis et alj|2002l) . due to its 
strong mid-IR absorption band centered at 9.6pm. On Earth, 
O3 is an indirect biomarker because its presence is due to the 
large O2 amount in the atmosphere, which itself is produced by 
the biosphere. Note that our model atmospheres do not contain 
O2 which would affect the temperature structure in our model 
mainly through its contribution to the heat capacity (Eq. [1} and 
the spectrum by a pressure broadening effect on spectral lines 



and potential collision-induced absorption, similar to N2, but the 
effect would be rather small. 

The following assumptions are used for the atmospheric 
composition: 



- CO2, O3 and N2 are assumed to follow isoprofiles. While 
CO2 and N2 are indeed isoprofiles (at least on Earth, Venus 
and Mars, where neither species significantly condenses), 
assuming an isoprofile for O3 is a simplifying assumption 
made in this work. Because O3 is formed by a 3-body reac- 
tion that requires both a pressure which is large enough and 
atomic oxygen produced by UV photolysis at low pressure 
levels, its abundance exhibits a very distinct maximum at 
mid-altitudes, as confirmed by many photochemical model 
studies. The terrestrial O3 profile can be approximated by a 
Gaussian profile without taking into acc ount interactive pho- 
tochemistry (e.g., Ivon Paris et al. 1 I2T5TT1) . Instead of one pa- 
rameter for the O3 isoprofile, this would introduce at least 
three parameters into the model, namely the location of the 
maximum, the concentration at the maximum and the full 
width at half maximum of the Gaussian which would com- 
plicate retrieval further (see also discussion of the effect on 
spectra in Sect. 12. 3t . 

- H2O and O3 are minor species, i.e. are only present at a level 
of several 10~ 2 or less, hence are neglected in the calcula- 
tion of the heat capacity or the mean molecular weight of the 
atmosphere. 

- Up to ptrop, H2O is calculated based on a fixed, constant 
relative humidity h, i.e. the partial pressure of water is ob- 
tained from Ph 2 o(T) = h ■ /? V ap,H 2 o(r), with the saturation 
vapor pressure according to local temperature. At p tmp , we 
assume a cold trap, i.e. for pressures lower than p tr0 p, H2O 
remains at the value calculated for ptrop- Using a constant 
relative humidity h in our model is also a simple assump- 
tion. On Earth, h is a function of altitude, longitude, latitude 
and highly vari able (standard globally-aver aging ID models 
often use, e.g.. lManabe & Wetheraldl[l967l which decreases 
monotonically with altitude). However, since the hydrolog- 
ical cycle on Earth is already difficult to model correctly, it 
will be much harder on exoplanets where no topographic or 
other constraints are available. Therefore, assuming a fixed, 
constant relative humidity in a first step is justified. 

- N2 is a filling gas, i.e. its concentration is adjusted such that 
the sum over all species is unity at the surface. 



Hence, the four atmospheric constituents are fully described 
by the temperature profile and three additional parameters (h, 
C0 2 and O3). 

Note that the possibility to constrain surface conditions and 
molecular absorption bands could be limited substantially by 
the pr esence of clouds (e.g., iTinetti et al.ll2.006l iKitzmann et alJ 
1201 ll) . In that sense, the cloud-free model used in this work il- 
lustrates upper limits on detectability. 

In order to calculate the heat capacity in eq. Q] we used the 
respective values for N2 and CO2, taken at 300 K. The values are 
(in units of cal K _1 moF 1 ) , c Pi n, = 6.953 and c p co^ = 9.215. 

The atmospheric altitude profile is then constructed by as- 
suming hydrostatic equilibrium, given the mean atmospheric 
weight, temperature-pressure profile and surface gravity. The 
surface gravity is calculated from the planetary mass and radius 
(mp and rp, respectively). At a given value of trip, rp is obtained 
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from a mass-radius relationship for rocky planets of ISotin et al.l 
(120071) (M e , R B Earth mass and radius, respectively). 



Temperature profiles 



Water profiles 



rp_ 

Re 



nip 
Mm 



(4) 



All of the 8 parameters (4 T/p, 3 chemical, 1 planetary) are 
uncorrelated in our model. In reality, this is not the case. For 
example, the surface temperature is of course related to sur- 
face pressure and atmospheric composition via the greenhouse 
effect. In turn, the chemical composition is influenced by at- 
mospheric temperature through equilibrium and photochemistry 
and temperature-dependent surface processes (e.g., outgassing, 
weathering). However, for the purpose of an illustration of a 
retrieval algorithm, a (self-) consistent atmospheric modeling 
would be computationally too expensive. Self-consistency is in- 
teresting when it allows to decrease the number of parameters. 
In the case of cool planetary atmospheres considered here, that 
are far from the thermodynamic equilibrium, self-consistency is 
out of reach: it implies, among other things, 3D time-dependent 
modeling with a realistic treatment of clouds and topography. 
Not only this cannot be used for retrieval but in practice it also 
increases the number of model parameters. 

2.2. Target scenario 

For our analysis, we used modem Earth as the reference case for 
an inhabited planet. We chose model parameters as to approx- 
imately reproduce the modern Earth temperature profile (see 
Table |2j. The relative humidity was set to /i=0.5 throughout the 
troposphere, C0 2 to 3.55 x 10 4 and 3 to 10" 7 . Note that an 3 
concentration of 3.75 x 10~ 7 would reproduce the atmospheric 
column amount of modern Earth. The planetary mass mp was 
fixed at 1 M e (hence, the model planet has a radius of 1 R @ , see 
Eq.H. 



Table 2. Target model parameters used in this work 



Ttoa [K] 

r surf [K] 

Pttop [bar] 
Pdry [bar] 
humidity 
CO2 vmr 
O3 vmr 
m P [M ffi ] 



270 

290 

0.3 

1.0 

0.5 

3.55 x 10- 4 

10- 7 

1.0 



Fig. |2 shows the temperature-pressure and water-pressure 
profiles in comparison with a modern Earth reference case. This 
reference case was obtained with a consistent ID radiative - 
convective climate m odel, coupled to ph o tochemistry (for a de - 
tailed description, see lRauer et al 1 1201 ll or lGrenfell et al.ll201lh . 

Despite the simplicity of the model, the general agreement 
between the target profile and the modern Earth reference is 
rather good. Also, the water profiles agree qualitatively well, 
even though the target profile is calculated with a constant rela- 
tive humidity. 

2.3. Planetary spectra 

The line-by-line code MIRART-SQuIRRL dSchreier & BottgeJ 
l200l was used to calculate high-resolution synthetic ra- 
diance spectra of the model atmospheres described above. 




target mode 
modern Earth 




200 220 240 260 280 300 
T [K] 



1 0~ 7 1 0~ B 1 0~ 5 1 0~ 4 1 0~ 3 1 0~ 2 1 0"' 10° 



Fig. 2. Temperature and water profiles of the target scenario, in 
comparison to modern Earth. 



Line parameters wer e taken from the Hitran 2008 database 
dRothman et al.ll2009h . The zenith angle was fixed at 38°. To 
obtain disk-integrated emission spectra, the radiance spectra 
were multiplied by tt. For the purpose of this work, c ollision- 
induced self-continua of CO? (I Words worth et al.ll2.Qldh and N2 
(iBorvsow & Frommhold|[T986i Laffertv et alj|1996l) have been 
included. 

The resulting high-resolution emission spectrum for our tar- 
get model of Table [2] is shown in Fig. [3] It is compared to a 
spectrum obtained using global mean modern Earth profiles (see 
Sect. 12.2b . The water features (around the 6.3 /urn band and in 
the rotation band longwards of around 16yum) do not differ by 
much, indicating that the choice of relative humidity is indeed 
not critical. Some noticeable differences, however, occur in the 
optically thin window region between 8 and 1 2 fim (due to dif- 
ferent surface temperatures, modern Earth 288 K, our target sce- 
nario 290 K and the different O3 profile) and around 7.7 fim (due 
to a strong methane absorption band, with methane not being in- 
cluded in our model atmospheres). Near the 15 /im CO2 band, 
some differences in the spectra are apparent which are due to the 
slightly different temperature structure in the stratosphere (see 
Fig.EJ. 

For modern Earth, the O3 band probes the upper troposphere 
and lower stratosphere, at pressures of about 10~ 2 to 10 _1 bar, 
where temperatures are increasing with altitude. Hence, an iso- 
profile which would reproduce the atmospheric column amount 
of O3 results in less emission from the O3 band since the contri- 
bution function would peak closer to the tropopause. This would 
imply that assuming an O3 isoprofile would underestimate the 
emission of the O3 band. However, many atmospheric model- 
ing studies have shown that, e.g., planets orbiting around M and 
K stars would not show a strong strato spheric inversion (e.g., 
ISegura et al1l2003l iGrenfell et alll2007l) . In these cases, an O3 
isoprofile would actually overestimate the O3 band. Therefore, 
since we do not include the influence of the central star in our 
model (see above), assuming an O3 isoprofile is justified for the 
sake of model simplicity. 



2.4. Observational model 

The measured quantity for spectral observations is the number 
of photons in a wavelength bin. To translate the calculated high- 
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Fig. 3. High-resolution emergent spectrum for the target model 
of Table [2] and a c omparison with modern Earth (taken from 
lGrenfelletai1l201lh 



resolution spectra from Sect. 12. 3 1 into a number of photons, we 
place the hypothetical target planet at a distance of 10 pc. We fur- 
thermore assume an ideal detector (i.e., throughput and quantum 
efficiency of unity) with a 10 m 2 collecting surface (i.e., about 
3.6 m circular telescope aperture) and an integration time of 1 
hour. Then, the spectra were binned to equidistant bins, corre- 
sponding to different spectral resolutions (/?=5-100) at 10/im 
(i.e., bin sizes of 2-0.1//m), over a spectral range of 5-20 pm, 
thus yielding a noise-free spectrum If. 

When comparing observed spectra and models to estimate 
atmospheric or planetary properties, it is important to consider 
the signal-to-noise (S/N) ratio per bin of the observations. We 
take the noise source Ns to be constant over wavelength. This is 
of course an approximation, since many individual, wavelength- 
dependent terms contribute to the overall noise (e.g., photon 
noise of both star and planet, thermal noise of the telescope and 
instrument, zodiacal light etc.). However, this approach yields 
a simple tool to calculate synthetic observations independent of 
telescope and instrument configurations. We define Ns as 



#288 (A)) 
fN 



const. 



(5) 



where Z?288 is a blackbody emitter of one Earth radius (R @ ) with 
a temperature of 288 K at a distance of lOpc, Aq- lOyum and 
an arbitrary noise factor. 

In order to construct noisy spectra 1^, we added a Gaussian 
noise of amplitude — yJWs to If : 



Im = hi + % ■ A t 



(6) 



where i is the wavelength bin and %, drawn randomly from the 
standard normal distribution JVo.i . The S/N ratio per bin / is then 



(S/N)i 



A n 



(7) 



In order to reproduce the EChO and Darwin S/N ratio spec- 
ifications at the respective spectral resolution (see Table [TJ, we 
adapted the noise factor fy in eq. [5] accordingly (e.g., for EChO, 
/ai=5-10~ 3 at R-\Q to produce a mean S/N per bin of 5-6). Fig. 
|4] shows examples of the thus generated synthetic spectra. 



Fig. 4. Example spectra of the target model (see Tableland Fig. 
for different mission designs as summarized in Table Q] 

2.5. Fit models 

When analyzing (noisy) observations and comparing them to a 
model, one standard quantity which is usually evaluated is the 
sum of the weighted squared residuals, the so-called^ 2 . The^f 2 
is defined as 



X 



Oi - Mix) 



(8) 



with x the vector of parameters (in our case, 8 parameters, see 
Sect. 12.11 ). Oi the noisy observation, M, the parameter-dependent 
model value and cr, the noise level (obtained from eq.|5j i.e. cr, = 
An — V^Vs" ) in a respective spectral bin ; and N is the number of 
data points. 

The^ 2 as defined above in eq.[8]is the log-likelihood in the 
case of Gaussian measurement errors, which we assume here. 
The ncr uncertainties on retrieved model parameters are then cal- 
culated by x 1 ~ Xmi„ ^ " 2 > e -g- a 3 cr uncertainty corresponds to 
a distance of A% 2 =9 from the minimum^ 2 value CtjL,). Finding 
a best-fit model corresponding to a certain minimum Xmi n > how- 
ever, does not guarantee a good fit corresponding to the global 
minimum. Two quantities are used generally to check the qual- 
ity of a fit, i.e. the p value and the reduced ^ 2 ed . The p value 
is the probability of generating a worse (i.e., larger) x 2 than the 
one calculated if the chosen parameters were the true ones. All 
parameter vectors x with^ 2 (jc)<,^ 2 would be considered a good 
fit, with^ 2 the^ 2 value corresponding to the chosen p. A stan- 
dard threshold is p=0.01, i.e. the false-alarm probability is 1 %. 
Based on the^ 2 , the^ 2 ed is calculated via 



2 , , X 2 (X) 



(9) 



where dj is the number of the degrees of freedom obtained from 
df = N — P (P number of parameters). In order to be considered 
a good fit, x 2 ed should be of the order of unity. 

In this work, we used two different approaches for com- 
paring modeled spectra with the synthetic observations. Firstly, 
to fit the observat ions, we used the IDL fitting routine MPFIT 
dMarkwardtl 120091) . which is an implementation of the nonlin- 
ear least-squares Levenberg-Marquardt (LM) algorithm dMorel 
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119781) . MPFIT uses the x" value to calculate best-fitting param- 
eters. However, since it is based on a Newton-type method, it 
might find a local minimum which could be quite far from the 
actual global minimum of the^ 2 . Note that Newton-type solvers 
are to some extent dependent on the initial parameter guess. 

Secondly, we calculate the x 2 f° r a large number of param- 
eter combinations to produce a x 2 ma P of the parameter space. 
With these maps, local minima in the^ 2 and potential degenera- 
cies between model parameter s could eas i ly be v i sualized. Such 
a method was used, e. g., by ILee et al. ( 20121). Cochran et~ai1 
(1201 lh . iFabrvckvetail (120121) . or lFord et all (1201 2b in addition 
to formal fitting procedures. Error bars or uncertainties on re- 
trieved parameters were then calculated from the x 2 maps. 

For a large number of parameters (e.g., 8 in the full model 
described in Sect. 12. Il l, producing x 2 maps is not practical, both 
in terms of calculation time and visualization. In such cases, the 
LM algorithm (or other algorithms, such as Monte-Carlo) is es- 
sential to estimate parameter values. 



3. Results 

3.1. Surface conditions 

In this section, we show the^ 2 maps in the r sul f-/?dry plane, i.e. 
the parameters relevant to surface habitability. The maps repre- 
sent slices of the parameter space with all other 6 parameters 
(planetary mass, atmospheric composition, upper atmospheric 
T/p profile) held constant at the values stated in Table [2] i.e. 
the true parameters used when generating the synthetic spec- 
tra. This means that uncertainties on surface conditions cannot 
be assigned based on such slicing, since for real observations, 
the model parameters are not known exactly a priori. However, 
these^ 2 maps illustrate the difficulty when trying to infer surface 
conditions from observations. The^ 2 values were calculated for 
20 synthetic observations, each with the specified S/N ratio per 
bin, but with a different realization of the noise. The median x 2 
of these 20 observations is plotted to illustrate the general shape 
of the x 2 ma P at the chosen S/N ratio. Note that if we were to 
combine these 20 observations (equivalent to adding transits), 
we would increase the S/N ratio by a factor of V20. 

Fig. [5] shows the^ 2 maps for the EChO (upper) and Darwin 
(lower ) low resolution con figuration, i.e. a spectral reso lution of 
#=10 dTinetti et al.ll2012l) and #=5 dLeger et al.lll996l) . respec- 
tively. The mean S/N ratio per bin is approximately 5-6 in the 
upper part and 12 in the lower part, i.e. close to the respec- 
tive S/N goals stated in Table Q] The different colors in the x 2 
maps represent increasing confidence levels, i.e. 1, 2, 3 and 5 <x 
confidence levels, corresponding to a distance of 1, 4, 9 and 25 
to the minimum y 2 . (see Sect. 12.5b . For EChO, at a 3 cr level 

^min v 1 ' ' 

(i.e., blue area in the x 2 maps), surface temperatures can be in- 
ferred to be higher than 280 K and lower than 300 K, indicating 
a habitable surface with relatively high confidence. However, for 
the extraordinary claim of the detection of (inhabited) habitable 
planets, one would need rather the 5 cr uncertainty level. This 
barely excludes surface temperatures below the freezing point 
of water, i.e. 273 K in the EChO case, while surface habitability 
can be inferred at the 5 cr level for the presented Darwin speci- 
fications. In contrast, surface pressure is less constrained in both 
cases (about a factor of 4), but still within reasonable values for 
habitability. 

Fig.|6lsho ws the y maps fo r the medium resolution D arwin 
cases (#=20. iLegeret al.lll996l and ICockell et al.ll2009al) . The 
mean S/N ratio per bin in these cases is 20 and 10, respectively 



>5 cr 
5 a 
3 a 

1 a 



>5 a 

5 a 
3 a 

2 a 
1 a 



Fig. 5. y 2 maps for EChO dTinetti et~ai1l2012l upper panel) and 
Darwin dLeger et al.lll996l R—5) specifications. Symbols show 
considered grid points. 



(see stated S/N goals in Table [T]). It is clearly seen that the rel- 
atively high spectral resolution together with the high projected 
S/N ratio per bin allows for a secure >5 cr characterization of 
surface conditions (but keep in mind the idealized setup of the 
X 2 maps in this section). 

In summary, in the idealized case presented here, the re- 
trieval of surface conditions is possible with a high S/N ratio 
(of the order of 5 and more). The 1 cr uncertainties are rela- 
tively small, but the optimistic case considered here underesti- 
mates retrieval uncertainties. Adding free parameters such as at- 
mospheric composition, stratospheric temperature structure and 
planetary mass to the fit will only increase the uncertainties of 
the retrieval. Furthermore, it should be kept in mind that the 
cloud-free nature of our model implies that retrieval possibili- 
ties for surface conditions might be over-estimated. 

3.2. Atmospheric composition 

The detection of atmospheric species is another major goal for 
exoplanet characterization. Habitability, in terms of surface con- 
ditions, does not necessarily imply an inhabited planet, hence the 
search for atmospheric biomarker species holds the greatest po- 
tential for actually finding extraterrestrial life. Therefore, in this 
section, we show the x 2 maps for the atmospheric composition 
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Fig. 6 . x 1 m aps for Darwin spec (R =20, upper: iLeger et all 
[19961 lower: ICockell et all l2009ab . Symbols show considered 
grid points. 



(i.e., h, CO2 and O3), again, as above, with all other parameters 
held constant at the true values. As before in Sect. 13. 11 the differ- 
ent colors in the x 2 maps represent increasing confidence levels 
and the median^- 2 of the 20 synthetic observations is plotted. 

Figure|7]shows the^ 2 maps for the EChO resolution {R-\Q) 
in the respective CO2-/1 (upper) and CO2-O3 (lower) plane. The 
S/N ratio is chosen to be equal to the EChO goal of 5 (see Table 
Q]). It is clearly seen that the presence of water can be securely in- 
ferred from the spectrum, as h-Q is excluded at high confidence. 
Also, CO2 can be detected clearly, since the minimum CO2 mix- 
ing ratio is about 10~ 5 . O3 is detected at the 1 cr level, but O3 
concentrations of 10~ 10 are compatible with the observations at 
the 2 cr level, hence the detection would be somewhat marginal. 
Note however that even 2 <x constraints would be very useful for 
the potential selection of candidates which would be subjected 
to further, more detailed follow-up observations. 

However, the goal of exoplanet space missions (such as 
EChO, Darwin, Finesse) is not only the detection of atmospheric 
species (i.e., strong lower limits), but also the characterization 
(i.e. a quantitative estimate of the concentration, meaning strong 
upper limits). As implied by Fig. [7] 5 cr uncertainties on CO2 
concentration cover about three orders of magnitude, and the hu- 
midity is only constrained to about a factor of 9 (with h-\ and 
h-Q excluded at high confidence). Therefore, it seems unlikely 



Fig. 7. x 1 contours for ECho specifications from Table [T] White 
symbols show considered grid points. 



that a S/N ratio of about 5 at spectral resolution R-IQ will allow 
for an accurate characterization of Earth-like atmospheres. 

In Fig. [8] we show the^ 2 maps with the same spectral reso- 
lution {R-\Q) but with a S/N ratio sufficient to limit the 5 cr un- 
certainties of the atmospheric composition to below one order of 
magnitude. The S/N ratio per bin (~28) is five times higher than 
the value stated in Table [1] but allows for the characterization of 
H2O, CO2 and O3 in this ideal model exercise (all parameters 
except composition fixed). An increase of S/N ratio from 5 to 28 
translates into a factor of about 30 in integration time if read-out 
noise is not the dominant noise source, because then, S/N scales 
with the square root of the integration time. 

In Figs. [9] and [TO] we show analogous x 1 maps for 
the Darwin low-resolution (R-5, see Table Q} and medium- 
resolution (R-20, see Table [TJ case. Since at low resolution, the 
aim is to constrain CO2, we only show the CO2-O3 plane in Fig 
[9] Results in Fig. [9] suggest that a S/N ratio of around 12 is not 
sufficent at R-5 to constrain CO2 better than about 2 orders of 
magnitude. To obtain one order of magnitude uncertainty at the 
5<x level, an increase in S/N ratio to more than 25 is needed 
(lower part of Fig . [9} . 

At medi um spectral r esolutio n of R=20, ICockell et al.l 
(l2009al) and ICockell et all (l2009bl) state that an S/N of 10 is 
aimed for to characterize the atmosphere with respect to CO2, 
H2O and O3. Hence, we show both the CO2-I1 (left column) and 
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Fig. S-x 2 contours for ECho resolution, with increased S/N ratio Fig. 9. x 2 contours for Darwin low resolution specification (up- 



of 28. White symbols show considered grid points. 



per) and needed S/N ratio (lower). White symbols show consid- 
ered grid points. 



the CO2-O3 (right column) plane in Fig. [10] As can be inferred 
from the upper panel (S/N ratio of 10) in Fig.fTOj except for CO2, 
atmospheric composition cannot be constrained with these spec- 
ifications. An increase to a S/N ratio of about 20 (lower panel) is 
needed to better constrain H2O and O3 at the 5 cr level (although 
O3 is already constrained at the 3 cr level for an S/N ratio of 10). 
Note that this S/N rati o at R—20 corresp onds to the stated O3 
characterization aim of lLeger et al.l(ll996l) in Table Q] 

An increase of about a factor of 2 in S/N ratio then translates 
into a factor 4 for Darwin in terms of exposure time necessary 
for a specific target. This significantly reduces the number of 
potential target stars. 

To further illustrate the difficulty to obtain firm constraints 
on O3 concentration, we show in Figs. fTTI and fT2l sample noisy 
observations together with spectra of model atmospheres con- 
taining various amounts of O3, but otherwise the 7 remaining 
parameters fixed at the "correct" values. Fig. [TT] shows a syn- 
thetic observation at R-\Q and a mean S/N ratio per bin of about 
6, i.e. the current EChO specifications. It is clearly seen that in- 
creasing the O3 concentration from 10" 10 to 10 6 changes the O3 
band (within 2 cr in the bin), however, the change in overall^ 2 is 
not significant, as expected from Fig [7] 

In Fig. |T2j the effect of O3 is shown with the same S/N ratio 
of about 6, but at a spectral resolution of R-IQQ. As expected, 



the higher resolution clearly improves upon the accuracy for O3 
retrieval. 

Therefore, our results imply that the current Darwin de- 
sign is probably overly optimistic with respect to the number 
of habitable terrestrial planets which could be characterized 
and searched for atmospheric biomarkers with emission spec- 
troscopy. Even putting aside the current paucity of targets, since 
no transiting HZ planet are known so far and only a few can- 
didate super-Earths in or near the HZ have been discovered (see 
Introduction), EChO is probably also unlikely to constrain atmo- 
spheric composition with the current S/N and R goals. 

3.3. Complete retrieval of planetary parameters 

In this section, we performed a non-linear least squares fit to 
one specific synthetic observation using MPFIT (see above, Sect. 
12.5b . All 8 model parameters were allowed to vary, as is the case 
for real observations, and we initialized the fit with a random 
guess of the parameters. Upper and lower boundaries for the fit 
parameters are stated in Table|4] The numerical values for the pa- 
rameter boundaries are mostly motivated by the fact we aim at 
investigating habitable, terrestrial planets, which, e.g., puts con- 
straints on the planetary mass and the surface temperature. In 
total, 100 randomly initialized fits were performed. 
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Table 3. Found best-fit parameters after 100 randomly initialized fits, with associated 3cr uncertainties. Spectral resolution and 
S/N ratio as proposed for EChO and Darwin, as well as increased S/N ratios. Parameters which could be fitted reasonably in red. 
Statistics of the fit results are also included (see Sect. 12.5b . 



parameter 


^=10, S/N=5 


R=\Q, S/N=28 


R=5, S/N=12 


R=5, S/N=26 


R=20, S/N=10 
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Fig. 10. x 1 contours for Darwin high resolution specifications 
(upper) and needed S/N ratios (lower). White symbols show con- 
sidered grid points. 



Fig. 11. Influence of O3 on synthetic observations (upper panel, 
noisy spectrum in black) and corresponding x 1 values (lower 
panel, diamonds). EChO resolution and S/N ratio. 



Table 4. Parameter range used for initial guess and retrieval iter- 
ation. 



parameter 


range 


Ttoa [K] 


200 - 350 


T^a [K] 


200 - 400 


p ttop [bar] 


0.01-0.75 


Pity [bar] 


0.5 - 25 


h 


0-1 


CO2 vmr 


- 0.99 


O3 vmr 


0-0.01 


m p [M e ] 


0.01 - 10 



Minimizing the x 2 determines the best-fitting, optimal pa- 
rameter vector x op t corresponding to xtm- The 100 fit results 
were used to determine jc opt . Uncertainties were then estimated 
based on the distance Ay 2 t-oxtm, ( see Sect. 12.51 ). 

Figs. [T3l and fT4l illustrate how this is done. They show the 
X 2 values as a function of the respective parameter values for all 
points evaluated by MPFIT during the fitting procedure (i.e., not 



only the best-fit values), a total of about 20,000 points. We chose 
the R-5, S/N=12 case, i.e. the low-resolution Darwin aim (see 
Table [T]). Indicated by horizontal lines are the Ay 2 =l, A;^ 2 =4, 
Ax 2 =9 and A^ 2 =25 thresholds, corresponding respectively to 
1, 2, 3 and 5cr uncertainties in the estimated parameter value. 
Also clearly seen are the local minima for T SUT f, CO2 and O3. 
Interestingly, the high O3 concentrations correspond to essen- 
tially zero CO2, indicating a strong correlation between both 
species at low spectral resolution. Furthermore, this also sug- 
gests that CO2 cannot be detected in this scenario, since no firm 
lower limits are found (see below, Table|3]l. 

As stated in Sect. 12.51 a best-fit model found by a fitting 
method not necessarily corresponds to a good fit. Therefore, 
Table [3] shows the quality of the fits as indicated by the p value 
and the;^ 2 ed value. Except for the R=5 observations where there 
are as many parameters as there are data points (d /=0, hence by 
definition^ 2 d =°°, Eq. [9jl, the^ 2 d indicates that good fits are in- 
deed obtained by MPFIT. Also, as stated in Sect.O MPFIT is 
dependent on initial conditions. This is indicated in Table [3] by 
the fact that less than 100% of the found best-fits are actually 
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Fig. 12. Influence of O3 on synthetic observations (upper panel, 
noisy spectrum in black) and corresponding y 1 values (lower 
panel, diamonds). /?=100 and S/N ratio of 6. 



Fig. 14. x 1 values as a function of parameter values: composition 
and mass, for R-5, S/N=12 (Darwin). Horizonal lines indicate 
1, 2, 3 and 5 <x uncertainties. 
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Fig. 13. x 2 values as a function of parameter values: temperature 
profile, for R=5, S/N=12 (Darwin). Horizonal lines indicate 1, 
2, 3 and 5 cr uncertainties. 

good fits, when adopting a p value of 0.01. However, generally, 
most estimates of jc op t actually are good fits, even though we 
used a Newton-type algorithm. 

In Table [3] we present the found parameters as well as the 
3<x uncertainties (i.e., Ay 2 =9, see Figs. [T3l and [T4li for several 
combinations of spectral resolution and S/N ratio. We use mis- 
sion specifications as indicated by Table[T]as well as simulations 
with increased S/N ratio. 

Table [3] shows that the temperature profile is largely uncon- 
strained by the observations. This is apparent from the found 
values of 7toa and p trop which are, within 3 cr, compatible with 
the boundaries imposed in Table|4] Only at R=20 with a S/N=20, 
fit results seem to indicate at least lower limits for 7toa- 

Limits on surface temperatures are tighter, with the fit re- 
sults in most cases excluding cold surfaces with high confidence. 
However, for the nominal EChO and Darwin R-5 cases, sur- 
face temperatures around 260 K are compatible with the data, 
which is inconsistent with surface habitability. Even at R-20 and 
S/N=10, surface temperatures above 273 K are only marginally 



inferred from the fit. Additionally, for the Darwin R-5 and the 
EChO nominal cases {R=5, S/N=12 and tf=10, S/N=5, respec- 
tively), upper limits for r sur f are consistent with surface temper- 
atures of the order of 400 K, wh ich is higher than the limit of ex - 
tremophilic life on Earth (e.g., iRofhschild & Mancinellill200ll) . 
Surface pressure is relatively well constrained to within a few 
bar at R-20 and, with an increased S/N of 28, also at R=10. 
For the Darwin R=5 cases, however, the entire pressure range of 
Table |4]is compatible with the noisy observations. In summary, 
surface habitability in terms of surface temperature could be in- 
ferred in most of the cases, but the temperature profile remains 
unconstrained. 

The planetary mass is generally well-fitted to within a factor 
of two or better to the true value, indicating that, combined with 
further constraints based on discovery data, it is probably not a 
critical parameter in our model. Of course, this is due to the fact 
that mass correlates with radius (see eq. |4), which then directly 
influences the emitting surface, hence the number of photons de- 
tected. Based on our assumed mass-radius relationship (see Eq. 
0), the uncertainties in planetary radius are of the order of 10- 
50 %, which is of the order of current uncertainties in transit 
surveys. 

With respect to atmospheric composition, Table [3] implies 
that the most difficult parameter for the retrieval is actually 
the humidity h. Most cases only allow for weak lower limits, 
indicating that a significant detection is missing. Also, giving 
firm upper and lower limits to characterize the humidity is only 
marginally possible for the R-20, S/N=20 case. 

The CO2 concentration cannot be constrained for the EChO 
nominal case, and even a detection is only marginal. At an in- 
creased S/N of 26, the EChO resolution allows for CO2 charac- 
terization within two orders of magnitude. For the Darwin R-5 
cases, only upper limits to CO2 could be provided, with lower 
limits not allowing for a conclusive detection. At R-20, in con- 
trast, clear upper limits are found. O3 cannot be characterized 
accurately with either R=5 or R= 10 nominal S/N ratios, but upon 
increasing the S/N ratio at the EChO resolution, relatively strong 
upper limits are found. A spectral resolution of R-20 allows for 
relatively firm upper O3 limits, however at S/N=10 a detection 
is rather marginal. 
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These results imply that even with increased S/N ratio (with 
respect to the current mission design, Table [T), both Darwin and 
EChO are unlikely to be able to provide accurate constraints on 
atmospheric composition and temperature structure, even though 
surface temperature and, to some extent, surface pressure are in 
principle accessible, confirming results from Sects. 13.1 l and [3~2l 

In Fig.[l5j we present scatter plots of the best-fit parameters 
at both the projected Darwin and EChO designs from Table Q] 
to investigate possible correlations between fit parameters. The 
most obvious correlations are between CO2, O3 and p& ly (upper 
panels). CO2 and O3 concentrations are generally lower at higher 
surface pressures. This is related to the depth of the 15 pm and 
9.6yum bands, respectively. The apparent correlation between 
CO2 and O3 is interesting in regards of the discussion on pos- 
sible false-positive detections of both species (see also discus- 
sion below). Correlations between CO2 and O3 with 7toa and 
Ptrop become apparent (middle panels), at least at better S/N ra- 
tio and higher spectral resolution (red crosses in Fig. Q3}. At 
higher values of p tmp , CO2 and O3 concentrations are systemat- 
ically lower. This is related to the parametrization of the model 
stratosphere. Both species are best probed in bands originating 
in the lower stratosphere. Since higher p lmp values mean warmer 
stratospheres (all other parameters being equal, see e.g. eq.[3]and 
Fig.Q]), higher p aof> therefore translates into lower concentrations 
of CO2 and O3. An analogous argument is valid for the appar- 
ent weak correlation between CO2 and O3 with 7Yoa- As stated 
above, there also is a weak correlation between planetary mass 
and r sur f. However, r sul f does not show any strong correlation 
with surface pressure or atmospheric composition (bottom pan- 
els). 



4. Discussion 

4.1. Atmospheric model 

In reality, a full atmospheric model which is able to reproduce 
observations of terrestrial, potentially habitable planets should 
contain many more parameters than the one presented here (see 

Sect. ITU- 
One important factor is the presence of clouds in the at- 
mosphere which have been omitted in the present work (see 
Sect. I2.1l i. Clouds could alter the atmospheric temperature pro- 
files and h ave a potentially large influence on surface tempera - 
tures (e.g jForget & Pierrehumbertll997llkitzmann et al 
Furthermore, they also a ffect the spectral app earance of a planet 
in emission spectra (e.g.. iTinetti et al.ll2006 ). They could mask 
molecular absorption bands and inhibit the probing of the lower 
atmosphere and surface (hence, the retrieval of surface condi- 
tions). On Earth, many different types of clouds exist in terms 
of, e.g., vertical and horizontal cloud extension, optical proper- 
ties or composition. Cloud cover is highly variable both spatially 
and temporarily. This adds many complications both for atmo- 
spheric modeling and spectral retrieval since usually, many sin- 
gle observations must be added to obtain reasonable S/N ratios. 

In addition, many more radiatively active species besides 
the ones considered in this work could have an impact on 
temperature structure and emission spectra. In particular, the 
biomarker species nitrous oxide could be detectable at least in 
high-resolu tion spectra given sp ecific planetary scenarios (e.g., 
early Earth. iGrenfell et al]|201 lb . SO2 and other sulphur species 
have been proposed in anoxic atmospheres or atmospheres 
with high volcanic activity (e.g . jKaltenegger & Sasselovll2010l 
iDomagal-Goldman et al. H2, given appropriate concen- 

trations of larger than »;10~ 3 , could also be detectable due 



to distinct collision-induced absorption (e.g.. lBorvsowl[T99Tl) . 
Additionally, the impact of the stellar radiation on the temper- 
ature profile should be taken into account since variations even 
for fixed atmospheri c composition are potentially large (e.g., 
iKitzmann etal]|2010l) . 

Furthermore, many photochemical studies showed that the 
influence of central star on atmospheric composition cannot be 
neglected. As an example, the biomarker-relevant gases methane 
and chloromethane could build up to large concentrations given 
a suitable stellar UV radiati on environment suc h as on planets 
around M dwarf stars (e.g., ISegura et al .1120051) . The CO con- 
centr ation was shown t o increase for planets orbiting M dwarfs 
(e.g. jRauer et al EoTTt) . However, using (photo-)chemical mod- 
els introduces many more additional parameters, e.g., number of 
important species for catalytic cycles, planetary surface fluxes, 
stellar and orbital parameters, geochemical cycles, etc., which 
would complicate the retrieval of planetary properties. Note, 
however, that constraining CO2 concentrations to within a few 
orders of magnitude might be indicative of the existence of 
a sili cate weathering cycle which controls CO2 (lAbbot et al.l 
|2012|) . Still, the potential importan ce of using c on sistent mod- 
eling t ools is demonstrated by, e.g.. ISelsisI (120001) or lRauer et al.l 
(1201 ll) . Coupling temperature structure, stellar radiation and at- 
mospheric photochemistry could cause biosignatures such as the 
O3 band to disappear in the emission spectra. 

Planetary and orbital properties such as rotation rate, obliq- 
uity and eccentricity induce diurnal and seasonal variability in 
the emission and reflectance spectra, which again poses a prob- 
lem for the retrieval, both because of additional free parameters 
and, in the case of thermal spectra, due to severe degeneracies 
between obliquity and th ermal inertia (e.g. jGaidos & Williams! 
12004 ICowanetal 1 l2012l) . Surface properties, e.g., oceans, con- 
tinents, vege tation, ice and snow cover, also influence planetary 
spectra (e.g. jHu et al Numerous studies have shown how 

to retrieve such information fr om spectrophotometry of Earth in 
the visible spectral range (e. g. JOaklev & Cashll2009l iFuiii et al.l 
l2010ilLivengood et al.l20llh . However, in a full planetary model 
such properties would increase the number of free parameters, 
hence render retrieval more complex and most likely prohibit 
conclusions on habitability. 

The bulk properties of a planet, e.g. its mass-radius rela- 
tionship, does not necessaril y follow the one o btained for rocky, 
terrestrial planets used here (ISotin et al. 2007). The masses and 
radii of super-Earths discovered so far span a large range (e.g., 



Leger et al.ll2009llCharbonneau et al.ll2009l.lLissauer et al.ll201 ll 



Batalha et al. 1201 lh . Assuming different bulk compositions for 
the p lanet then leads to v ery different mass-radius relationships 
(e.g- lFortnev et al.ll2007h . 

Regarding the parametric nature of our model, note that self- 
consistent modeling could provide a better constraint on CO2 
and, indirectly, on O3, since detecting O3 has a lot to do with con- 
straining CO2 (see discussion below). Therefore, the retrieval re- 
sults presented here might be conservative. Future studies could 
use real observations of Earth or Mars, or results from consistent 
modeling studies, as target scenarios and investigate whether 
accurate and correct retrieval is possible in such more realistic 
cases. 

4.2. Observational constraints 

In order to reduce the number of fit parameters and break exist- 
ing degeneracies in emission spectra, observational constraints 
and additional information must be used when interpreting the 
spectra. Many important planetary properties could be deter- 
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Fig. 15. Scatter plots of best-fit parameters from 100 randomly initialized fits. Black: EChO case, R-10, S/N=6, red: Darwin 
medium-resolution case, R=20, S/N=10. 



mined beforehand, such as mass and radius. Measurements of 
the planetary mass are ava ilable through a nu mber of meth- 
ods, such as astrometry (e.g., Malbet et al.l201 ll) . radial velocity, 
transi t timing variations in multip le systems (e.g. jLissauer et al 



201 1 ) or dynamical models (e.g. jBeust et al 



( e.g..lLiss auer 
. 2008, Mayor 



et al 



2009). The planetary radius could be estimated for transiting 



planets. 

In addition, for transiting planets it is possible to com- 
bine emission spectroscopy with transmission spectroscopy. 
Transmission spectroscopy is, to first order, sensitive to at- 
mospheric scale height, hence mean atmospheric weight and 
planetary gravity, and only weakly dependent on the verti- 
cal temperature structure. Thus, degeneracies between temper- 
ature and atmospheric co mposition could be broken (see, e.g., 
iBenneke & Seagerll2012l) . Due to the large spectral range cov- 
ered by EChO, such transmission observations could be per- 
formed with the same mission as the emission spectroscopy, in 
contrast to Darwin which would need to rely on other obser- 
vation facilities. Another complementary method, especially on 
non-transiting planets , could be to use p hase curves and the vari- 



ation spectrum (e.g.. Ha rrington et 



se pnase curves ana tne van- 
aLll2006l ISeager & Demingl 



l2009LISelsiseTal1l201 1, Maurin et al. T l2012h . a technique which 



could be used by the Darwin misson. 

Broadband optical photometry (over a spectral range of 
about 0.2-2 //m) could be used to measure the reflected light 
Ir of the planet, which would impose constraints on the Bond 



albedo As and the planetary radius (/« ~ Ab ■ r 2 p ). Since Ag 
depends strongly on atmospheric composition and surface pres- 
sure through absorption of stellar radiation and Rayleigh scatter- 
ing and Eq. prelates mass and radius, such measurements could 
in principle provide useful additional constraints. However, the 
contrast in reflected light between star and planet is of the or- 
der of 10 8 -10 10 for terrestrial, habitable-zone planets. Hence, 
observations would be very challenging and, consequently, con- 
straints actually may be rather weak. 

The orbital period of the planet is determined from the 
discovery data rather precisely. The orbital distance is deter- 
mined from the orbital period based on stellar mass (and/or 
radius, for transiting planets). Combining orbital distance with 
the stellar energy distribution obtaine d either from spectrosco py 
or stellar model atmo spheres (e.g., iBuser & Kuruczl U992I or 
lHauschildt et al.|[l999l) . it is possible to compute the stellar en- 
ergy input on the planet. Then, assuming global energy conser- 
vation for the planet, it would be possible to eliminate model 
scenarios which would be prohibited, thus possibly limiting the 
parameter space. 



4.3. Real observations 

In practice, many of the aforementioned potential limits on pa- 
rameters are rather difficult to obtain. For example, the detec- 
tion and verification of a planetary signal is very challenging for 
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low-mass planets (see, e.g 
201 ll lHatzes et al l 1201 ll 



. iFerraz-Mello et al.ll201ll 



f or the CoRoT-7 system, 



Pont et al 



Vogt et al 



20 1 01 lTuomill20 1 ll lGregorvll20 1 ll for the GL581 system). Also, 



uncertainties on stellar parameters, usually derived from stel- 
lar models, are the main problem for derived planetary param- 
eters such as ra dius and mass, as is for example the case for 
GJ 12 14b (e.g.. ICarter etaDl201 ll ISada et al.ll2010l iBerta et alJ 
1201 lb where the estimated planetary radius changes by 15 %, 
depending on the adopted stellar model. 

Other practical concerns are instrumental constraints such 
as available filters, wavelength coverage and spectral resolution 
of the telescopes used for the observations or large time spans 
between observations. Constraints obtained with the same mis- 
sion (or instrument, even) are usually preferrable to constraints 
obtained with different observational setups due to systematic 
noise. 



4.4. Biomarker and bioindicators 

In this work, we considered ozone as a biomarker, since on Earth 
it is mainly produced from oxygen photochemistry, and the oxy- 
gen is provided by the biosphere. 

However, interpreting the presence of ozone as a bioindi- 
cator is not necessarily straightforward. In other words, a 
2 cr detection of ozone is not equivalent to a 2 cr detection 
of life. Many studies have shown that ozone can be pro- 



duced abiotically (e.g., iSelsis etaD 120021 ISegura et al.l 120071 
iDomagal-Goldman & Meadowsl2010l) . Additionally, it has been 
shown that ozone persists in the atmospheres of Earth-like plan- 
ets even at drastically reduced oxygen concentrations and that 
ozone concentrations d epend crucially on t he UV radiation fiel d 
of the central star (e.g. JSegura et al.ll2003UGrenfell et al.ll2007h . 
Very thin ozone layers have indeed be en detected both in the 
atmospheres of Venus and Mars (e.g., iLebonnois et al.l 120061, 
iMontmessin et al.l201 ll) . although concentrations are far too low 
to being detectable on exoplanets. 

The absence of ozone inferred from observations, however, 
does not imply that ozone is not present in the atmosphere, 
as a masking of absorption features is certainly possi ble, as 



sh own by, e.g..lSchindler & Kastind (120001) . ISelsis et all (120021) 



or lvon Paris et al.l (1201 ll) 



The most promising way of identifying biospheric signatures 
in atmospheric spectra seems therefore the simultaneous pres- 
ence of species s uch as water , ozon e, oxygen and methan e, as 
proposed by, e.g JSagan et al.l(ll993l) or lSelsis et al.l (120021) . 



5. Conclusions 

In this work, we investigated the possibility of the retrieval of 
planetary properties (such as atmospheric composition, temper- 
ature structure, radius) from emission spectra. We used a para- 
metric atmosphere model to create arbitrary atmospheres of hy- 
pothetical habitable planets. The atmospheric profiles were then 
used to calculate synthetic, noisy observations. Based on these 
observations, we tried to fit the model parameters and investi- 
gate potential degeneracies. 

We have used two design concepts (in terms of S/N ratio and 
spectral resolution) of exoplanet space missions. Results imply 
that surface conditions can be characterized relatively well (to 
within ~10K at 3o") with S/N ratios between 10-30, depending 
on spectral resolution. This would allow for an assessment of the 
potential habitability of the target planet. However, even with 
high S/N ratios of the order of 20 or higher, emission spectra 



did not allow for accurate determination of atmospheric compo- 
sition. The biosignature of O3 as well as the atmospheric H2O 
content are not well constrained with current mission designs. 
The temperature structure could not be retrieved. 

Since achieving such high S/N ratios (or even higher values) 
is rather difficult in the foreseeable future, single-visit emission 
spectroscopy alone is most likely not capable of characterizing 
the atmospheres of potentially habitable planets. In order to ob- 
tain higher accuracy on retrieved parameters of habitable plan- 
ets, emission spectroscopy must be combined with other tech- 
niques, e.g. transmission spectroscopy, phase curves or photom- 
etry in the visible spectral range. An investigation of such pos- 
sibilities, and whether this would allow for atmospheric charac- 
terization, is subject of future studies. 
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